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Abstract. A possible definition of the specific heat of open quantum systems is based on the reduced 
partition function of the system. For a free damped quantum particle, it has been found that under certain 
conditions, this specific heat can become negative at low temperatures. In contrast to the conventional 
approaches focusing on the system degree of freedom, here we concentrate on the changes induced in the 
environment when the system is coupled to it. Our analysis is carried out for an Ohmic environment 
consisting of harmonic oscillators and allows to identify the mechanism by which the specific heat becomes 
negative. Furthermore, the formal condition for the occurrence of a negative specific heat is given a physical 
interpretation in terms of the total mass of bath oscillators and the system mass. 
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1 Introduction 

The standard definition of the thermodynamic equilibrium 
quantities in a canonical ensemble relies on the assump- 
tion that the coupling between the system of interest and 
its thermal environment is arbitrarily weak. In the real 
world, this approximation never strictly holds and impor- 
tant deviations from the usual picture can be expected in 
particular for systems on the nanoscale. Beyond the study 
of effects of finite coupling to the environment on the ther- 
modynamic properties of system, the more fundamental 
question arises of how to properly define quantities like 
an internal energy, a specific heat or an entropy in such a 
situation [T]. The last five years or so have seen consider- 
able activities addressing this issue and related ones [2 ,3, 

In the quantum regime, the thermodynamic equilib- 
rium properties not only depend on the strength of the 
coupling to the environment. Even more importantly, it 
turns out that the very definition of these quantities is 
no longer unique [1] . As an example, we consider here the 
specific heat which for a damped quantum system can be 
defined in at least two different ways. 

The first approach defines the internal energy of the 
system as expectation value of the system Hamiltonian 
taken with respect to the thermal ensemble of the cou- 
pled system. By the latter we mean here and in the fol- 
lowing the coupled complex of system (S) and heat bath 
(B). Quantities related to the coupled system will be in- 
dicated by a subscript "S-l-B" . Based on the definition of 
the internal energy just introduced, one obtains the asso- 
ciated specific heat by taking the derivative with respect 
to temperature. Obviously, in the absence of a coupling 
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between system and bath, the specific heat defined in this 
way reduces to the standard specific heat of the system. 

An alternative approach is based on the reduced par- 
tition function of the system defined as pUll^E^ 



Z = 



--S-I-B 



(1) 



In the absence of a coupling between the system and its en- 
vironment, Z equals the partition function of the system. 
Even though we thus can consider 2^ to be a quantity as- 
sociated with the system, we omit a subscript "S" in order 
to emphasize the fact that Z also takes into account the 
coupling between system and bath. Here, we will employ 
the ratio ([I]) only in the context of equilibrium thermody- 
namics, but recently it has also been used in the discus- 
sion of fluctuation theorems for open quantum systems in 
nonequilibrium situations [71I51IT7]. 

The partition function ([T]) can be used to define ther- 
modynamic equilibrium quantities by means of the stan- 
dard relations valid in the absence of a coupling between 
the system and its environment. We thus obtain a second 
expression for the specific heat reading 



C = kB0 



2d ' HZ) 



(2) 



where /3 — l/fceT. In the absence of a coupling between 
system and bath, this definition of a specific heat also 
reduces to the usual definition and thus agrees with our 
first definition in the limit of vanishing coupling. Again we 
omit a subscript "S" to underline the dependence of the 
specific heat on the coupling, but we do so also in view of 
the interpretation ([s]) given below which emphasizes the 
role of the heat bath. 
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It should be noted that, in general, the two specific 
heats just defined differ for finite coupling to the environ- 
ment already in the leading high-temperature corrections 
to the classical specific heat. For a bilinear coupling be- 
tween system and bath, the classical specific heats still 
agree [HIS] while this need no longer be the case for an- 
harmonic couplings [S]. 

The difference between the two approaches is most 
spectacular in situations where the specific heat ^ be- 
comes negative in contrast to the positive specific heat 
based on the expectation value of the system Hamiltonian 
[l][3] . In the following, we exclusively consider the specific 
heat obtained from the reduced partition function so that 
the symbol C will always refer to ([2|. 

Although the reduced partition function ([ij appears 
to be a rather formal starting point, the quantities de- 
rived from it have a clear physical meaning. Due to the 
fact that thermodynamic equilibrium properties depend 
on the logarithm of the partition function, these quanti- 
ties in the presence of a finite coupling to a heat bath 
possess a natural interpretation in terms of the difference 
between the quantity evaluated for the coupled system 
and the quantity evaluated for the bath alone [Sll23j. For 
the specific heat ([2]), we thus have 



C = Cs+B — Cb ■ 



(3) 



sisting of several harmonic oscillators is coupled to a finite 
environment |19| . 

In the treatment of dissipative systems, the system 
degree of freedom usually is in the focus of interest. On 
the other hand, the difference ([3| of specific heats, even 
though conceptually associated with the system degree of 
freedom, provides the motivation to take a different point 
of view, namely the one of the environment. The question 
answered by the specific heat ^ really is: How does the 
specific heat of the environment change when a system 
degree of freedom is coupled to it? 

In the following we shall address this question by study- 
ing the change in the spectral density of the environment 
when a system degree of freedom is coupled to it. In Sec- 
tion [2] we briefly introduce the model for the damped free 
quantum particle which will allow us to obtain the change 
in the spectral density of the environment in Section [3] 
From these results, we derive in Section |4] the thermo- 
dynamic properties of the system in the sense explained 
above. We will show that these results can indeed be ex- 
pressed in terms of properties of the system degree of free- 
dom. In Section [5] we will give a physical interpretation of 
the condition under which the specific heat (|3 1 of the free 
damped quantum particle becomes negative. Finally, we 
present our conclusions in Section [6j 



The two quantities involved in this difference, Cs+b and 
Cb, constitute measurable quantities. Even though one 
may invoke a superbath imposing the temperature on the 
bath and on the coupled system, the coupling to the su- 
perbath may be assumed to be negligible. Therefore, we 
are certainly allowed to employ the standard thermody- 
namic definition for the quantities on the right-hand side 
of ([3]). C thus is not a specific heat in the proper sense 
but it is the natural replacement for a specific heat in 
a situation where the system of interest is coupled with 
non-negligible strength to an environment. The approach 
discussed here was recently employed in the analysis of 
thermal data obtained for a metal phosphate compound 
[23] . For the sake of simplicity, we will occasionally refer 
to C as specific heat of the system. In doing so, we should 
however keep in mind the preceding discussion. 

With the above interpretation of the specific heat ([2|, 
negative values should not give rise to concerns with re- 
spect to thermodynamic instabilities. In fact, there is no 
reason why the difference of two positive quantities should 
be positive. The situation considered here should therefore 
be distinguished from the negative specific heat appear- 
ing e.g. in self-gravitating systems within a microcanonical 
description pS p g ]. !??]. 

Apart from the free damped quantum particle which 
we focus on in this paper, negative entropies and/or spe- 
cific heats have been discussed e.g. in the context of the 
Casimir effect |28ll291 l30 l [3"i p 32 p33ll341 l35;. Kondo systems 
[551157] . XY spin chains [H], two- level fluctuators [5] and 
energy transport in proteins [38j . While for a single har- 
monic oscillator coupled to an Ohmic bath ^j|3j or to a 
single bath oscillator [B] the specific heat was found to be 
positive, this is not necessarily the case if a system con- 



2 Free particle coupled to harmonic oscillators 

As our model for the study of the appearance of a nega- 
tive specific heat we consider a particle of mass M which 
is bilinearly coupled to a set of harmonic oscillators consti- 
tuting the environment. The corresponding Hamiltonian 
is given by 



H = 



2M 



E 



pI 

2m„ 



(4) 



The sum describes the environmental oscillators of mass 
m„ and frequency w„, the bilinear coupling and a poten- 
tial renormalization. The latter ensures that the effective 
equation of motion of the system position Q corresponds 
indeed to the Langevin equation of a free damped particle 
taking the form 



MQ + M / ds7(i - s)Q{s) = i{t) . 



(5) 



Below, we will relate the damping kernel ^{t) to the prop- 
erties of the environment. The properties of the noise term 
f (<) are irrelevant for our discussion so that we refer the 
reader to the literature for details [20,21,,,22.;39]. The cou- 
pling constant of the bilinear term in Xn and Q in Q 
has been expressed in terms of the masses and frequencies 
of the environmental oscillators which is possible without 
loss of generality [40,41 . The Hamiltonian is manifestly 
invariant under spatial translations of all degrees of free- 
dom, confirming once more that we are treating a free 
damped particle. 



Gert-Ludwig Ingold: Thermodynamic anomaly of the free damped quantum particle: the bath perspective 



3 



In our discussion, two spectral densities play a key role 
and it is important to carefully distinguish them. We refer 
to the first quantity as spectral density of eigenmodes p. 
This quantity is defined in terms of the eigenfrequencies 
of the quadratic Hamiltonian Q. For the bath alone, we 
have the spectral density of eigenmodes 



/5b(w) 



n=l 



(6) 



Correspondingly, we define a spectral density of eigen- 
modes ps+b('^) in which the frequencies cj„ appearing in 
(|6| are replaced by the eigenfrequencies of the coupled 
system described by Q. The difference ps+B — Pb wiU 
form the basis for the derivation of various thermodynamic 
quantities in Section [4] 

In addition, in order to characterize the properties of 
the heat bath and its coupling to the system, a quantity 
commonly referred to as spectral density of bath oscilla- 
tors \M 



J{lu) 



IS 

n=l 



(7) 



is introduced. In contrast to the spectral density of eigen- 
modes Q in the absence of coupling, the spectral den- 
sity of bath oscillators ([t]) not only depends on the fre- 
quencies of the bath oscillators but also on their masses. 
For our choice of coupling constants (cf. discussion below 
([5|) this is tantamount to saying that ([7| depends also on 
the coupling constants between system degree of freedom 
and bath oscillators. J{u)) contains all information about 
the heat bath required to describe the properties of the 
damped quantum system. 

Although it is not obvious, it will turn out that even 
within our bath-centered approach all quantities of inter- 
est will eventually be expressible in terms of the spectral 
density of bath oscillators ([t]) . This already indicates that 
even though we take the point of view of the environment, 
we will learn something about the properties of the sys- 
tem. 

In our results, it will be natural to express the depen- 
dence on J{lu) in terms of the Laplace transform of the 
damping kernel ^{t) appearing in ([5]). The two quantities 
are related by means of [20.,21j.22j 



7(z) 



2 



z 



Q W -I- 



(8) 



For the purpose of our discussion, we will restrict our- 
selves to so-called Ohmic damping which implies that the 
spectral density of bath oscillators is continuous and linear 
in the frequency at least for small frequencies. At high fre- 
quencies, the spectral density of bath oscillators will typ- 
ically be suppressed with respect to this linear behavior. 
To be specific, in cases where such a cutoff is of relevance, 
we will employ the so-called Drude model with 



(9) 



Here, the low-frequency behavior is characterized by the 
damping constant 7 and uj-q is the cutoff frequency. The 
corresponding Laplace transform of the damping kernel 
reads 

7d(^) = — • (10) 



Because the damping constant 7 combines with the in- 
verse temperature /3 to form a dimcnsionless quantity hj3^ 
and thus merely sets the temperature scale, the cutoff fre- 
quency wd will be an important parameter determining 
the properties of the damped free particle. 



3 Change in the spectral density of bath 
oscillators 

In our environment-centered approach, we start out with 
a set of uncoupled environmental oscillators with frequen- 
cies Un- In addition, we have one system degree of free- 
dom corresponding to an undamped free particle. We now 
couple the environmental oscillators to the free particle 
as prescribed by the Hamiltonian Q. As discussed in the 
previous section, the system and environment together are 
translationally invariant. 

In the spirit of Ullersma's analysis of the damped har- 
monic oscillator [43j we determine the eigenmode spec- 
trum of the free particle coupled to its environment. Diag- 
onalizing the Hamiltonian, we recover a zero-mode which 
replaces the zero-mode corresponding to the uncoupled 
free particle. We thus concentrate on the non-zero frequen- 
cies which due to the coupling are shifted with respect to 
the original environmental oscillator frequencies. The new 
frequencies are obtained as solutions i7 of the equation 



C30 9 



122 



= M . 



(11) 



For the following intermediate steps it is convenient to 
explicitly consider an equidistant set of discrete frequen- 
cies UJn — nA with n = 1,2,... for the environmental 
oscillators. Later, we will take the limit of vanishing fre- 
quency spacing A to recover a continuous spectral density 
of environmental oscillators. In view of 0, we have the 
relation 

dc. J{u:)f{u) ^\Y. ^niolfinA) (12) 

71=1 

which holds for any function / for which the integral and 
the sum exist. If, on the other hand, we replace the integral 
on the left-hand side by a Riemann sum with step width 
A, we can express the masses m„ of the environmental 
oscillators in terms of the continuous spectral density of 
bath oscillators J{uj) according to 



2 J{nA) 
TT {nAy 



A. 



The eigenfrequency condition (11) thus becomes 



E 



J{nA) 



-M . 



(13) 



(14) 
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In order to treat frequencies fl close to an unperturbed en- 
vironmental frequency nA correctly, we rewrite this con- 
dition as 



1 



n=l ^ ' 



J{nA) Jif2)\ _ IT 



(15) 



i7 / 



M 



Performing the first sum exactly and replacing the second 
sum by an integral, we finally obtain 



cot 



7ri7 



A 



with 



(16) 



(17) 



The coupling of the environment to the system degree 
of freedom modifies the original spacing A between adja- 
cent environmental eigenfrequencies yielding a new spac- 
ing A + e. The correction e depends on the spectral bath 
density via the function g defined in ( 17 ) and can be shown 
to be of order Exploiting the latter fact and making 
use of the addition theorem of the cotangent, one deter- 
mines e from (16). For the change in the spectral density 



of eigenmodes, one thus finds as a central result 



Ps+b(^)-Pb(^?) = 



1 

A + e{[2) 



1 

A 



1 g'm 

TT 1 + .g(i7)2 



(18) 



Here, the prime denotes a derivative with respect to the 
argument. 

Further insight into the change of the spectral density 
of eigenmodes can be obtained by specifying the properties 
of the bath. We choose a Drude model with the spectral 
density of environmental oscillators given by Q and a 
Laplace transform of the damping kernel as specified in 
(lol). From (Il8| one obtains 



Ps+B — Pb 



where 



UJ2 



l?2 +, 



f2^ 



Wl.2 = Cli' ± iuj" 



~2 



l± Wl 



1^2 



47 



(19) 



(20) 



and is the Drude frequency providing the high-fre- 
quency cutoff. These three frequencies are the eigenfre- 
quencies associated with the deterministic version of the 
Langevin equation ([s]) for a Drude damping kernel [5]. 

The change in the spectral density of eigenmodes is 
shown as thick solid line in Fig. [l] for various values of 
ujYill- The dotted lines represent the Lorentzian contri- 
butions according to (19). For the sake of clarity, they 



are depicted also in the grey regions of negative frequen- 
cies which are irrelevant for the spectral density of eigen- 
modes. As long as the cutoff frequency is sufficiently high. 



(a) 
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Fig. 1. The change in the spectral density of eigenmodes 
Ps+B — Pb induced by the coupling to the system degree of 
freedom is depicted as thick solid line for Drude damping with 
(a) cjd/7 = 5, (b) u)-a /7 = 1, and (c) uioM ~ O-l- The contribu- 



tions of the three Lorentzians according to ( 19 1 are represented 
as dotted lines. Only the white regions of positive frequencies 
are of relevance for the spectral density of eigenmodes. 



wd > 47, ( |19[ ) describes the sum of three Lorentzians cen- 
tered at zero frequency. An example for this situation is 
given in Fig. [l^. Coupling the system degree of freedom to 
the environment then leads to an increase of the spectral 
density at zero frequency by {uj-q — 7) /ttjloy) . For smaller 
cutoff frequencies, if wd < 47 like in the cases depicted in 
Figs. [Th and c, the two frequencie s ([20| become complex 
and the first two Lorentzians in ( |19[ ) are now centered 
at nonzero frequencies. As a consequence, the suppress- 
ing effect of the third Lorentzian on the spectral density 
of eigenmodes becomes relevant. For wd < 7 one finds in 
fact a suppression of the spectral density at zero frequency 
as is shown in Fig.[lj; for wd /7 = 0.1. It is this suppression 
which leads to the negative specific heat as we will discuss 
in the following section. 
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Before doing so, we would like to emphasize that al- 
though the change in the spectral density of eigenmodes 
obtained in this section is mainly related to properties 
of the environment, the results are determined by prop- 
erties of the damped system. In the general result ( 18 1 



the environment appears only through the spectral den- 
sity of bath oscillators Q which is sufficient to provide 
a reduced description of the damped system. In the spe- 
cific case of Drude damping, the three eigenfrequencies of 
the damped system turn out to be sufficient for the com- 
plete description of the change in the spectral density of 
eigenfrequencies . 

Despite being closely related to the properties of the 
damped system, the change of the spectral density of eigen- 
modes discussed here should not be confused with the den- 
sity of states which can be obtained from the reduced par- 
tition function ([!]) by its inverse Laplace transform [44] . 
Because the density of states defined in such a way is lin- 
early related to the reduced partition function, it cannot 
be interpreted as a difference of a property of the coupled 
system on the one hand and the heat bath on the other 
hand 



4 Thermodynamic properties of the system 

The change of the specific heat ^ due to the coupling of 
the system degree of freedom to the environment can be 
obtained from the change ( [l8| of the spectral density of 
eigenfrequencies as 



C = I da;(ps+B(w) - pB(w))Cho(t^) 
/o 



where 



Cho(w) 



hf3uj 



2smh{hl3u}/2) 



(21) 



(22) 



is the specific heat of a harmonic oscillator with frequency 



ui. Eq. 21 implies a negative specific heat C at low temper- 
atures if the spectral density of eigenmodes is suppressed 
for small frequencies due to the coupling of the system 
degree of freedom to the bath. 

Assuming Ohmic damping, i. e. J(w) ~ uj at low fre- 
quencies, we can obtain the specific heat at low temper- 
atures by inserting th e ch ange in the spectral density of 
eigenmodes (18) into (plj). Together with (17 1 one finds 



C 



TTl + f (0) kBT 



(23) 



3 7(0) h 

in agreement with the findings of Ref . [3l . For 

7(0)<-l, (24) 

we thus obtain a negative specific heat at low tempera- 
tures. A physical interpretation of this result will be given 
in Section [5] below. 



More explicit results can be derived for the case of 
Drude damping. There, the integral in (21) can be evalu- 
ated analytically and one finds 



C 
fee 



V 2tt 



V 27r 



27r 



(25) 

where ij^' {z) denotes the digamma function. This result 
is in agreement with the expression obtained by proceed- 
ing according to the point of view of the system [3] and 
starting with the reduced partition function ([T]). The two 
approaches therefore are equivalent, but the approach pre- 
sented here gives additional insight thr ough the change in 
the spectral density of eigenvalues (18). 



Instead of obtaining the internal energy and the spe- 
cific heat from the reduced partition function, we can use 
our result (25) for the specific heat to obtain the other 
two quantities. In the limit of vanishing coupling to the 
heat bath, the internal energy U is related to the specific 
heat C by means of 



C 



~d/3 



(26) 



In view of ([s]) we thus obtain from ( 25 ) by means of an 
integration the difference of internal energies induced by 
the system-bath coupling 



U = Us+B - Ub 



V 27r j 
~2^^ I 27r / 



~2^^ V 27r 
1 

2p- 



(27) 



This result is only determined up to a constant of inte- 
gration which, by comparison with the known result [3], 
turns out to vanish. In particular, at zero temperature, we 
thus find 



In 



WD 



tVjJ2 



In 



WD 
Ul2 



(28) 



Furthermore, by means of the relation between the in- 
ternal energy U and the partition function Z 



U 



d/3 



HZ) 



(29) 



one reproduces the correct temperature dependence of the 
ratio of the partition functions of system and bath on the 
one hand and bath alone on the other hand .3i 



1 



r I 



27r y V 27r 
2tt / 



r I 



(30) 
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This result leaves a prefactor undetermined which for the 
thermodynamical equilibrium quantities is irrelevant. The 
temperature dependence, however, agrees with the result 
obtained by means of other techniques, e.g. the path inte- 
gral approach. 



5 The missing mass of bath oscillators 

In the previous section, we have found that a negative 
specific heat occurs provided the condition (24 1 is satis- 



fied. This rather formal condition can be given a physical 
meaning. In order to avoid an infrared divergence in the 
limit of vanishing argument z, we first express the Laplace 
transform of the damping kernel ([s]) as 

ii^) = m + - T'i^i^-m)-,^- (31) 



Muj 



With (24) and (31) the condition for a negative specific 



heat then becomes 



M7(0)a; - J{uj) 



-/ dw ^^^^5 ^>M. 

TT Jo uj-^ 



(32) 



Observing that according to (It]) the total mass of the 
bath oscillators can be obtained from the spectral density 
of bath oscillators J(a;) as [iUllil] 



ri=l 



2 J(c.) 
nin = - / do; — — 



the condition ( 32 ) can be expressed as 



AMb > M . 



(33) 



(34) 



Here, AMb is defined by the left-hand side of ([32j) and 
refers to the total mass of oscillators which are missing in 
the actual bath with respect to a strictly Ohmic reference 
bath where J(a;) — M^{Q)ijj for all frequencies. We note 
that while due to the infrared divergence mentioned above 
the total mass is infinite for both baths, the difference 
AMb is finite. 

We thus arrive at the following interpretation. Cou- 
pling a system degree of freedom of zero frequency to the 
bath oscillators tends to shift the frequencies of the latter 
to larger values. However, for a strictly Ohmic environ- 
ment, i.e. in the absence of a finite cutoff in the spectral 
density of bath oscillators, the ensemble of bath oscillators 
resists a reduction of the spectral density of eigenmodes 
at low frequencies. A suppression only becomes possible 
if one eliminates a number of bath oscillator with a total 
mass at least as large as the mass associated with the sys- 
tem degree of freedom. The bath then is no longer able 
to resist the "pressure" of the system degree of freedom 
pushing the eigenmodes to higher frequencies. 



6 Conclusions 

In order to achieve a better understanding of the unusual 
thermodynamic properties of a free Brownian quantum 



particle, we have employed a somewhat uncommon ap- 
proach to the analysis of dissipative quantum systems by 
taking the perspective of the environment. This approach 
is, however, quite natural if thermodynamic quantities are 
defined in terms of a reduced partition function of the 
system because they actually refer to the change of these 
quantities when the system degree of freedom is coupled 
to the heat bath. 

We have analyzed the modification of the bath spec- 
trum induced by the coupling to the system degree of free- 
dom. Interestingly, the spectral density of bath oscillators 
is sufficient to describe the change in the bath spectrum. 
As a consequence the latter depends only on properties of 
the damped system as has been exemplified by means of 
a Drude-type damping. Starting from the change in the 
bath spectrum, expressions for the specific heat, the in- 
ternal energy and the reduced partition function obtained 
previously from a system-based approach have been re- 
produced. 

The low-temperature behavior of the specific heat of 
the free damped particle is determined by the shift of the 
low-frequency environmental oscillators induced by the 
coupling to the system degree of freedom. If the ratio 
of cutoff frequency and damping strength is sufficiently 
small, the system degree of freedom succeeds in suppress- 
ing the spectral density of eigenmodes at low frequencies 
by shifting the bath modes to higher frequencies. As a re- 
sult, at low temperatures the specific heat of the heat bath 
is lowered if the system degree of freedom is attached. In 
contrast, for larger cutoff frequencies, the high-frequency 
oscillators of the environment act against the tendency of 
the system to shift the bath modes to higher frequencies. 
Then, the specific heat remains positive for all frequencies. 

The condition to be satisfied by the environment to 
allow for a negative specific heat at low temperatures 
has been shown to have a physical interpretation. The 
anomaly in the specific heat appears if the mass miss- 
ing in the environment with respect to a strictly Ohmic 
reference bath exceeds the mass associated with the sys- 
tem degree of freedom. Then, the bath oscillators are not 
strong enough to resist the zero-frequency degree of free- 
dom which, when coupled to the environment, tends to 
increase their frequencies. Otherwise the spectrum of bath 
modes is sufficiently stiff to prevent the spectral density 
of eigenmodes from being reduced. 

The author is grateful to Michele Campisi, Peter Hanggi, Astrid 
Lambrecht, Serge Reynaud, Peter Talkner, and Juan Diego 
Urbina for useful discussions. This work has been supported 
by the DAAD through the PROCOPE program. 
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